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PREFACE 


This  is  the  first  report  in  a  series  of  studies  conducted  to  produce  a 
realistic  operations  research  theory  of  attrition  and  reliability  of  land 
supply  carriers.  This  research  is  a  part  of  the  Generalized  Transpor¬ 
tation  Analogue  (GTA)  effort.  The  purpose  of  the  "Truck-Transport 
Model"  study  is  to  test  conclusively  the  viability  of  (and,  if  viable,  to 
exploit)  the  hypothesis  that  general  terrain  properties,  enemy  action, 
personnel  behavior,  and  limited  command  decision  determine  the  most 
probable  performance  of  a  land  supply  system  and  that  this  performance 
can  be  adequately  described  with  transfer  theory. 

Chapter  1  states  the  basic  concepts  of  transfer  theory,  develops  the  sub¬ 
sidiary  terms,  and  derives  the  formulation  relating  to  these  terms.  Not 
explicitly  stated  is  the  assumption  that  the  processes  which  such  a  model 
describes  are  not  strictly  "mechanistic",  but  that  the  "on-going  events" 
are  affected  by  the  ensemble  of  past  occurrences. 

Chapter  2  explains  the  subsidiary  terms  employed  in  Chapter  1  and  re¬ 
lates  them  to  the  systems  research  problem  under  study.  Equations 
2.  6.  3  express  the  final  form  of  the  theory  and  are  amenable  to  the 
following  interpretation: 

They  determine  the  number  of  loaded  trucks  in  transit  per  unit  area 
at  a  given  place  at  a  given  time  which,  in  the  course  of  carrying  out 
a  mission,  are  being  subjected  to  being  put  out  of  operation,  deflected, 
or  forced  ter  change  their  speed  by  terrain,  enemy  action,  or  decision, 
barring  historical  accident. 

In  the  succeeding  chapters  of  the  report,  this  theory  will  be  tested  and 
applied  against  theoretical,  practical,  and  historical  evidence. 
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CHAPTER  1  -  THE  TRANSFER  EQUATIONS 


INTRODUCTION 


Today,  military  mobility  is  becoming  identified  with  a  major  problem 
area:  logistics  of  motor  transport.  A  further  simplifying  assumption 
frequently  made  is  that  this  problem  can  be  measured  by  measuring  the 
ability  of  a  vehicle  to  traverse  natural  terrain.  This  paper  contends  that 
the  validity  of  this  implicit  assumption  has  not  been  demonstrated;  rather, 
the  complete  system  should  be  analyzed,  with  simplifying  assumptions  made 
only  if  the  model  complexity  is  beyond  the  capacity  of  modern  mathematics 
and  computer  technology. 

This  paper  will  present  the  derivation  of  the  theoretical  model.  The  im¬ 
portance  of  the  parameters  will  be  the  subject  of  a  computer  analysis  after 
the  model  has  been  historically  verified.  The  companion  paper  "The 
Terrain  Scattering  Probability"  is  prepared  separately  so  that  the  require¬ 
ment  for  a  measure  of  performance  for  specific  units  may  be  met  but  will 
enter  this  model  as  the  scattering  parameter. 


GENERAL  NATURE  OF  THE  PROBLEM* 

It  is  desired  to  know  the  transfer  intensity  and  the  time-space  density  of 
units  of  supply  under  large-scale  operations  involving  varying  battle  con¬ 
ditions  and  terrain  features.  We  assume  that  general  coefficient  and  proba¬ 
bility  functions  that  take  these  factors  into  account  can  be  realistically  de¬ 
fined.  Here,  we  will  investigate  the  intensity  and  density  of  loaded  carrier 
units.  In  general,  such  an  ensemble  can  not  be  adequately  described  in 
ordinal-/  space,  but  requires  a  six-dimensional  place  space  p  (kj  ,  x2  ,  x3 , 
vx  1  >  vxz  >  vx 3  ) .  This  problem  requires  only  a  four-dimensional  place 
space,  which  will  be  modified  to  suit  the  problem. 


DEFINITIONS 

— >  — > 

Consider  a  place  space p  with  coordinates  p  (R,  f2,  t)  , 

— ) 

Variables:  Position  =  R 

Direction  =  0  ■  (||  Q  ||  =  1  ). 

Dependent 

Variable:  Iy  =  intensity. 

*This  development  parallels  Chandrasekhar,  Reference  3,  and  follows 
Anselone, Reference  6.  (References,  are  given  in  the  Bibliography.) 
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1. 2. 1  I  (R,  )  =  units  at  R  traveling  in  the  direction  ft 

per  unit  area  _L  to  ft 

per  unit  solid  angle  about  ft 

per  second 

where  the  units  are  loaded  trucks. 


) 


For  a  given  R  and  ft,  construct  an  element  of  area  do  _L  to  ft. 

With  each  point  of  da  as  a  vertex,  construct  a  solid  angle  dw  about  ft. 
The  resultant  figure  is  called  a  pencil. 


F igure  1 .  Pencil. 


— >  — > 


From  the  definition  of  Iv  (R,  ft  },  we  obtain 
— ^  — >  — > 

1. 2.  2  Iv  (R,  ft  )  do-  du  dt  =  units  at  R  traveling  through  points  of  da 
with  directions  in  dui  about  ft. 

I  *“>  j 

Let  da  be  an  area  not  necessarily  -L  to  12, 


The  radiation  passing  through  da  also  passes  through  da1.  Thus, 

1.2.  3  IV(R,  £2  )  cos  9  da  dw  dt  =  unittf' aj^  R  traveling  through  points 
of  da  with  direction  in  du  Jtbout  ft. 


It  is  convenient  to  define 

— ^  ^  — >  — )  — )  — y  — )  — > 

1.  2.4  Iy  (R,  0  )  =  Iv  (R,  )  fl  ,  which  is  a  vector  at  R  in  the  ft 

direction  having  a  length  equal  to  the  intensity  in  that  direction. 
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The  net  vector  flux  it  Fv  (R)  is  defined 

— )  — )  .  — >  — )  — >  — >  — )  — ) 

1.2.5  ttFv  (R)  =  Ju  Iv  •  nQ  n0  dw  =  Iy  Q  .  n0  n0  dw. 

The  net  transfer  of  units  across  a  unit  area  in  any  arbitrary  direction 
N  at  R  (  ||N  II  =  1)  is  given  by 


1.2.6 

ttFv  • 

N  =  fu  Iv  ■ 

— ) 

n0 

— ■>  — )  — )  — )  — >  — > 

n0  •  N  dw  =  /w  I  ft  •  n0  n0  •  N  dw 

— ) 

L 

and 

1  1  TT  Fy  ||  = 

I  v  cos  0  dw  • 

where  0  is  the  angle  between  the  chosen  direction  N  and  the  variable 
(of  integration)  direction. 

Consider  the  units  associated  with  Iv  (R,ft)  for  a  particular  R  and  ft. 
As  the  units  proceed  in  the  ft  direction,  they  may  be  absorbed  or 
scattered  in  a  direction  or  velocity  other  than  the  direction  and  velocity 
of  inte rest. 

Thus,  the  intensity  is  diminished  as  it  travels  through  the  medium. 


Figure  3. 


— >  — >  — >  — )  — ) 

Let  Iv  (R  +  sft,  ft)  denote  the  number  of  units  of  Iv  (R,  ft)  that  remain 

— ) 

in  the  ray  after  traveling  a  distance  s  in  ft  direction. 

The  difference 

dlv  =  Iv  (R  +  ds  ft,  ft)  -  Iv  (R,  ft)  <  0 


is  generally  proportional  to  ds,  to  Iy  ,  and  to  the  density  p  of  the  medium. 
Thus,  write 


1. 2.  7  dlv  =  -  ky  p  Iv  ds 

where  ky  is  the  "Mass  Absorption  Coefficient". 

— ?  — ^ 

Consider  the  loss  from  Iv  (R,  ft)  due  to  absorption  in  an  element  of 
"mass  ", 

dM  ~  p  da  ds  (see  Figure  3). 
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Recalling  the  definition  of  I v (R,  A)  and  equation  1.2.2,  ,  the  losses  of 
units  in  dM  are  given  by 

1.2.8  Iy  da  •  ky  p  ds  du>  dt  =  kv  Iv  dM  dw  dt  =  loss 

of  units  of  velocity  v,  in  direction  Si  in  dM 

— ^  ^ 

Suppose  that  a  unit  traveling  in  tlm  A_direction  is  absorbed  at  R. 

We  define  the  Phase  Function,  p(A,  A'),  such  that 
— )  — > 

1.2.9  ^ ^  -  =  the  probability  per  unit  solid  angle  about  A'  that 

— ) 

the  unit  will  be  Scattered  into  the  A  direction. 


Thus, 

— ^  — N 

1. 2.  iO  J  ^  f1  >  ky  Iv  (R,  A')  dt  dM  dw  doi' 

->  -> 

=  units  scattered  from  direction  A'  into  direction  A  . 


Integrating  over  all  directions  A  1  of  the  incident  ray, 

1.2.  11  jv  /u,  p  (A,  A')  Iy  (R,  A1)  dco  1 

— ^ 

=  units  scattered  from  all  directions  A1  into  the  particular 
direction  A 

— ) 

per  unit  solid  angle  about  A 

per  second 

per  unit  mass  at  R  . 

{  g  ) 

jv  is  the  emission  coefficient  due  to  scattering.  The  general 
emission  coefficient  is 


1.2.12  jv  =  jv  (R,  A)  =  units  emittedMn  direction  A 
per  unit  solid  angle  about  A 
per  second 
per  unit  mass  at  R. 

The  ratio 


r  r  — *  — * 

1.2.13  j  ky  =JV  =/v  (R,  A)  is  the  source  function. 
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DERIVATION  OF  THE  ELEMENTARY  EQUATION  OF  TRANSFER 

_ ) 

Consider  an  element  of  mass  dM  emitting  in  the  ft  direction. 


dm  =  pda  ds 


Figure  4. 

From  the  definition  of  emission  coefficient,  we  have 

— > 

1.3.1  j  v  dM  =  j  v  pd  cds  -  units  emitted  by  dM  in  the  ft  direction 

-•> 

per  unit  solid  angle  about  ft  ,  per  second. 

Dividing  by  dcr  and  letting  dff->0,  we  obtain 

— )  —)  — )  — "4 

1.3.2  j  ypds  =  units  emitted  along  the  line  (R,  R  +  ds  ft  )  in  the  ft 
direction 

— ) 

per  unit  solid  an^le  about  ft 
per  unit  area  1ft 
per  second. 

By  definition, 

jv  p  ds  -  intensity  of  direction  ft  emitted  along  the  segment 
(R,  R  +  ds  ft  ). 


Taking  into  account  absorption  and  emission,  the  net  change  in  the  intensity 
is 

1. 3.  3  d  Iy  (R,  ft  ). 

Thus,  the  transfer  equation  is 

1.3.4  ±l^LA±.  =  Iy  (R,  ft  )  -  lv  (R,  ft  ) 
or 

1.3.5  J—  ft  •  VIv  (R,ft  )  Iv  (R,  ft  )  -  jy(R,  ft). 

^  d  — * 

Remember  that  ^  and  V  are  operators  for  the  variable  R  and  arbitrary 

but  fixed  ft  . 
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THE  DENSITY  OF  UNITS 


The  density  Uv  of  units  of  the  radiation  at  velocity  v  at  any  given  point 
is  the  number  of  units  per  unit  volume  that  are  in  transit  in  the  immediate 
neighborhood  of  the  point  considered. 

Following  Chandrasekhar  (Reference  :3,  Chapter  1,  Section  2.  3)  to  find 
this  density  at  a  point  P  ,  construct  around  P  an  infinitesimal  volume  V 
with  a  convex  surfacea  .  Around  V  construct  another  convex  surface  2 
such  that  the  linear  dimensions  of  2  are  large  compared  with  those  of  a 
and  such  that  the  volume  inclosed  by  2  is  nevertheless  small  enough  such 
that  the  intensity  in  any  given  direction  can  be  considered  as  the  same 
for  all  points  inside  2. 


Now  the  radiation  traversing  the  volume  V  must  have  crossed  some 
element  of  the  surface  2  .  Let  d2  be  such  an  element.  Further  let  0 
and  £  denote  the  angles  which  the  normal  to  d2  and  to  an  element  da 
of  a  makes  with  the  lines  joining  the  two  elements.  F rom  equation  1.2.3, 
the  units  streaming  across  d2  that  also  flow  across  da  are 


Let  l  be  the  length  traversed  by  the  pencil  of  radiation  through  the 
volume  element  V  .  The  radiation  1.4.1  incident  on  da  per  second  will 
traverse  the  element  in  a  time  ~ ,  where  v  is  the  speed  of  the  units.  The 
contribution  to  the  total  amount  of  radiating  units  in  course  of  transit 
through  v  by  the  pencil  is 


i.4.2  i  cosi_cos_0dodZ  l=  1  i  dvdoo 

V  rZ  V  V  v 


where  dw  =  dZ  is  the  solid  angle  subtended  by  dZ  at  P  and 

dV  =  da  cos  £  is  the  volume  intercepted  in  Y  by  the  pencil  of  radiation. 
Thus,  the  total  number  of  units  in  course  of  transit '  through  x  due  to- 
streaming  in  from  all’directions  can  be  obtained  by  integrating 
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SOURCE 


1.4.  2  over  all  V  and  w  .  Hence, 
1.4.3  i  /  dVf  du  Iv  =  ^/Iv  dw. 


Thus, 

1.4.4  IV  =  i  Jlv  du>. 

Let  Jv  =  /  Iv  dw  be  the  average  intensity.  Then 


1.4.5  IV  = 


4ir  T 
v  Jv 


NOTES  ON  THE  SPECIFIC  TRANSPORT  EQUATION 


EVENT  AND  INTERACTION  CROSS  SECTIONS 
(See  Figure  6) 


Absorption  (losses) 

Annihilation  by  enemy  action 

Speed  transfer  (enemy  action) 
Speed  transfer  (terrain) 
Scattering  (enemy  action) 
Scattering  (terrain) 

Incidental  losses 
Gains 
Sources 

Speed  transfer  (enemy  action) 

Speed  Transfer  (terrain) 

Scattering  (enemy  action) 
Scattering  (terrain) 

Incidental  gains 


kavi  Pavi 

k  av2  P av2 

k  bvi  P  avi 

kbv2  Pavz 

k  cvi  Pbvi 

k  CV2  Pbv2 

kdvi  Pavi 

kdv2  Pav2 

k  evi  Pbvi 

k  ev2  Pev 2 

S,  S2 

kbvz  Pavz  Pa  (v2  It,  '  ) 

kbvi  Pavi  Pa  (vi  -»v2  »  S2 ,  fi 1  ) 
kcv2  Pbv2  Pb  (v2  ^vi»  ) 

k  cvi  Pbvi  Pb  (vi  ~ )V2  >  fijfi1) 
kdvi  Pavi  Pcvi  (fi>fi')  kdv2  Pav2  Pcv2  (fl,fi') 

k evi  Pbvi  Pdvi  .)  k ev2  Pev 2  Pdv2  ) 

G1  G2 


where  the  coefficients  and  probabilities  are  defined  in  the  obvious  manner 
with  appropriate  dimensions. 
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(A)  Absorption: 


and  similarly  define 

^V2  ds  -  -  (  |  kav2  +  k  bv2  +  kcjv2  |  pav2  +  k  CV2  pbv2  + 

k  ev2  Pev  2  I^V2  ds . 


( B)  Emission  . 

j<JV]  (  v2  — >  V1  )  _  -  JGj'  Pa  (v2  — >  V1  >  ^ 1  )  f  V2  (R>  ^  )  dw'  > 

4  t r 

Kb )  f  — >  — )  — )  -4 

jflvi  (  v2  ->  V!  )  =  — - 1  V  pt,  (v2  _>v,  ,  12,  12'  )  Iv2  (R,  12 1  )  du>'  , 

r  4  TT 

K.^-y-  j  r  — >  -*>  — 4  — > 

jyvi  =  '4“" .  V  Pcvl  (  0.  O' )  IVI  (R,  Q'  )  du',  and 

j 6 vi  =-^-  L'  Pdvl  (H,  S'  )  IV1  (R,  fi'  )  du>\ 

4  TT 

Thus, 

1.5.1  d  IV1  (R,  12)  -  -  Kyl  IV1  ds  +  (Jqvi  \v2  vi  }  P av2  ^ 

J (3 vi  {v2  vi  J  Pbv2  d  j-yyj  Pavi  +  Jfivi  Pbvi  )  ds  ,  + 
Combining  similar  terms, 

1.5.2  d  IV]  (R,  12  )  =  -  KV1  IVJ  ds  +  j  v,  ( v2  v,  )  ds  +  j  VJ  ds 


where 

Jvi  (v2_*  Vj  )  =  (jay,  {v2->  v,  }  Pav2  +  jpvl  {v2  ->  V,  }  pbv2  ) 

and 

Jvi  _  0  yvi  Pavi  J6vi  P  bvi  )■ 

*In  the  case  of  v2  ,  terrain  scattering  density  pbvl  becomes  pcv2  . 
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Finally 


1.  5.  3 


dlvi  (R,R) 


ds 


~  -  KV1  ly]  (R  ,  R  )  4  jy]  (^2— *V|)  +  jy| 


For  convenience,  write: 

1  y-  — >  — >  — >  — >  ^ 

J  vi  =  Trr  f  u j1  |kdvi  Pa  vi  Pcvi  ( i2 »  )  +  k  evl  Pbvi  Pdvi  ^  1  )| 


— >  — > 


Iv,  (R,  n')dJ 


*V1 


j  vi  (v2  ->  vi  )  “  Tff  L'  fk  bv2  P av2  Pa  (v2  -*  vi >  ^ ^ 1 )  +  k  CV2  PbV2 
Pb  (V2  ->  V1  •  )}  I V2  (R>  ^')  dw'' 


Or 

J  vi  -  4tT  /w1  Pvi  (  f2 ,  R  1  )  Ivi  ( R,  fl  1  )  do 1  and 
j  VI  (v2  ->  Vj  )  =  jj  pvl  (v2  _>  v,  ,  R,  R!  )  Iv2  (R,  R'  )  dw'  . 

The  appropriate  systems  of  equations  now  take  the  form 
dl^  (R.fl) 


1.5.4 


1 


ds 


+  Kv[  IV1  (R,  R)  -  /w i  Pvi  (S2>  £2 1  )  IV1  (R,  R'  )dco' 
+  4^  f^1  Pvi  (v2  _>  Vj  ,  R,  R!  )  Iy2  (R,  R1  )  do! 


-4  — ) 

d  I  -  ( R  f  £2 )  — >  — )  i  — >  — )  — >  — > 

1.5.5  —  +  Kv2  I  v2  (R,  R)  -  /w'  Pv2  ( R .  R 1  )  V2  (R.  R')du' 


ds 


1 


+  4^  Pv2  (vi  ->  v2  >  ^2 ,  R 1  )  I  vi  (Ri  R  ')  dco1. 


These  equations  are  of  the  operational  form. 


1.5.6 

Dlvi 

=  kvi 

I'Vl 

+  AV1 

^  V2 

rsj 
> 
p — i 

Q 

n 

< 

tv. 

IV2 

+  AV2 

I VI 

Or 

Dx  = 

?? 

< 

X 

+  Avz  y 

Dy  = 

k  v2  y 

+  Av2  X 

• 
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So  that 


1.5.7 


AVi 

y  = 

Dx 

-  kVi 

X 

> 

< 

rj 

X  = 

Dy 

“  k  v2 

y 

This  is  the  general  model  without  arbitrary  source  or  Iobs  terms, 

THE  TIME-DEPENDENT  TRANSFER  EQUATIONS 


-»  -> 


(A)  In  the  simple  transfer  model,  if  Iv  (R,S2,  t)  is  time-dependent,  the 
transfer  equation  takes  the  form 

,  ,  .  1  dI(R,Q,t)  1  d Iv(R.  8,0  A>  M  ,  T,ff  -A  ti 

1-6-1  k^77  — irt —  i^Tp  — "Si -iv(R,n,t)  +  v(R*«.t) 

— )  — ) 

(B)  In  the  two-speed  transfer  problem,  if  Iy(R,  fl,t)  is  time-dependent, 
the  transfer  equations  take  the  form 

1.6.  2  —  d  --vl  (R-’  =-  dIvl  -  Xvj  I vi  (R.fl.t) 

vi  dt  ds 


+  Pvi  )  1V.  (R.O!,  t)  du' 

+  4Vf„'  PV1  (Vj  ->v,  ,  SJ,  n1  )  IM  (R,  !!’  ,t)  do,' 

i  div2  (R.n.t)  _  diV2  (R.n.t) 


■1.6.3 


~sr 


-  K  vi  1  vi  (R»  8  »  0 


+  L'  Pvz  0  dw' 

+  iV  b'  pvz  (V! ->  v2,S2,fi' )Iv,  (&»"&.  0  du>', 

where  arbitrary  source  and  loss  terms  are  deleted. 
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CHAPTER  2 


OPERATIONAL  SIGNIFICANCE  OF  THE  TRANSFER  MODEL 
INTRODUCTION 

Before  proceeding  either  to  a  theoretical  investigation  of  the  system  of 
equations  1.5.4  and  1.5.5  or  1.6.2  and  1.6,3  derived  in  the  last  chapter 
or  to  studying  them  in  a  specific  geometry,  we  will  relate  these  systems 
of  integrodifferential  equations  to  the  systems  research  problem  at  hand 
through  a  statement  of  the  significance  and  dimensions  of  the  coefficients 
and  probabilities  involved. 

Consider  a  large-scale  theater  of  operation  involving  regions  of  varying 
terrain  conditions  and  various  intensities  of  enemy  aclion.  Suppose  that 
it  is  desired  to  be  able  to  predict  or  maintain  a  supply  system  (function) 
in  this  area.  To  do  this  it  is  assumed,  barring  local  anomalies  and 
historical  accident,  that  by  and  large  the  gross  determination  of  a  supply 
function  is  governed  by  general  terrain  properties,  enemy  action,  personnel 
behavior,  and  command  decision.  In  this  paper,  we  will  consider  only 
supplies  delivered  by  land  carriers  such  as  trucks  and  will  define  our  unit 
of  supply  as  one  loaded  truck.  We  assume  that  this  unit  travels  at  a  speed 
v  and  is  destination  oriented.  In  the  course  of  a  "trip",  it  may  be  put  out 
of  operation  (annihilated),  deflected  (scattered),  or  forced  to  change  its 
speed  (transferred)  by  terrain,  enemy  action,  or  decision.  Taking  such 
effects  into  account,  the  problem  is  to  construct  an  applicable  model  that 
determines  a  probable  supply  function. 

We  have  accepted  the  fact  that  a  natural  dependent  variable  or  concept  is 
the  intensity  (defined  in  Equation  1.2.  1)  of  units  at  a  point  in  space  traveling 
in  a  given  direction.  A  realistic  formulation  of  absorption  and  scattering 
requires  this  angular  dependent  intensity  and  requires  an  integral  formu¬ 
lation  as,  for  instance,  the  inte  -ity  in  a  given  direction  at  a  point  is 
influenced  by  units  being  scatt  into  this  direction  from  all  other  di¬ 
rections  at  this  point.  We  are  thus  led  in  a  plausible  manner  to  a  transfer 
equation,  and,  because  of  its  generality,  adopt  transfer  theory,  which  we 
will  attempt  to  make  applicable. 

In  the  specific  model  considered,  it  is  further  assumed  that  cross-country 
mobility,  vt,  substantially  differs  from  on-road  mobility,  v2  ;  for  this 
reason,  two  speed  groups  are  considered.  This  leads  to  a  system  of  two 
dependent  transfer  equations,  one  governing  the  intensity  for  on-road 
speed  and  the  other  governing  the  intensity  for  off-road  speed  where  the 
average  speeds  are  taken  to  be  representative  of  the  groups. 

We  now  have  to  relate  the  coefficients  and  probabilities  occurring  in  the 
transfer  equations  to  operationally  meaningful  quantities  and  derive  re¬ 
lations  expressing  these  quantities.  The  coefficients  and  probabilities 
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will  be  given  for  a  two-spacial  dimensional  model.  The  preceding  model 
and  theory  have  been  presented  in  three-dimensional  form.  There  are 
two  virtues  in  this  approach.  The  first  is  an  emphasis  on  generality  and 
the  second  is  to  distinguish  clearly  the  nature  of  physical  entities  in  similar 
applications  and  the  operational  entities  that  will  be  presented  below.  For 
instance,  there  is  no  conceptual  difficulty  in  defining  our  two-dimensional 
density.  The  two-dimensional  equations  to  be  considered  are  similar  to 
those  already  derived  except  for  the  occurrence  of  a  factor  of  2  ir  in  place 
of  4  it  (see  equation  1.2.9)  with  the  obvious  changes  in  dimensions. 

GENERAL  SIGNIFICANCE  OF  THE  DENSITY  p  AND  MASS  ABSORPTION 
COEFFICIENT  k. 

Strictly,  the  density  p  -- 

’  7  p  "  A~*0  A 

where  A  is  an  area  and  M  is  the  mass  associated  with  this  area.  As 
we  are  interested  in  a  stochastic  model,  average  values  and  probabilities 
apply,  so  define 

_M_  lim  M 
Pz  ~A  Z  A->0  A 

where  M  is  the  total  mass  associated  with  a  fixed  area  A. 

— )  — )  — }  — ) 

Recall  that  I  (R,  fi  )  dadw  dt  =  units  at  R  traveling  in  direction  fl. 

and  klv  (R,  A)  dcadtdM  “  loss  of  units  at  R  traveling  in 

direction  A. 

Thus,  take  kpds  to  be  dimensionless. 

Now  p  =  mass/unit  area,  where  mass  is  defined  in  terms  of  intensity 
of  enemy  action,  terrain  properties,  etc. 

k  is  thus  seen  to  be  a  probability  that  in  a  unit  area  an  absorption  will 
take  place 

per  unit  mas s 
per  unit  length. 

k  can  be  called  the  mass  absorption  coefficient. 
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OPERATIONAL  DEFINITION  OF  THE  DENSITIES  AND  MASS 
ABSORPTION  COEFFICIENT  ~  . . ~ 


We  now  consider  separately  the  different  absorption  terms  given  in 
the  table  (page  10 ). 

Annihilation  by  Enemy  Action 

Define 

p3iVl  =  occurrences  unit  area,  where  an  occurrence  could  be  the 
detonation  of  a  bomb.  As  discussed, 

^  j  _  total  number  of  expected  events  (occurrences)  in  a  region 

^avi  “  area  of  the  region 

Then  kavl  is  the  probability  per  unit  occurrence  that  a  unit  will  be 
annihilated  in  terms  of  length  J_  to  ft  . 

aav2  and  kav2  are  similarly  defined. 

Speed  Transfer  by  Enemy  Action 

ky3V1  and  kbv2  are  just  the  probabilities  per  unit  occurrence  that  a 
speed  transfer  will  take  place. 

Speed  Transfer  by  Terrain 


The  nature  of  terrain  is  fundamentally  different  m  the  cases  of  on-road 
and  off-road. 

Off-Road:  From  terrain  features  such  as  slope  .and  soil  conditions, 
ago/no-go  function  can  be  evaluated  for  an  area  which  will  determine 
whether  a  specific  vehicle  type  will  or  will  not  stall  in  this  area.  If  we 
consider  a  region  R  (see  Figure  7)  and  impose  a  grid  on  it  giving  go/no-go 
values  for  different  cells, 


no-go 


go 


R 

Figure  7. 
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Then  k  becomes  the  probability  per  unit  no-go  that  a  unit  will  be 
transferred  given  in  terms  of  length _  to  d  . 

On  -Road  :  Consider  a  road  net  in  a  given  region.  Impose  upon  it 
a  rectangular  grid  (see  Figure  8)  such  that  the  map  in  each  cell  is  con¬ 
nected  and  such  that  each  crossroad  falls  within  a  cell.  Now,  in  each 
cell  of  this  grid,  construct  a  linear  approximation  to  the  road  map  such 
that  the  maximum  distance  from  the  road  to  the  linear  segment  is  a  pre 
assigned  value,  saye  (see  Figure  9),  and  such  that  the  segment  is  tange 


to  the  road  at  its  contact  points.  Now  consider  the  following  two 
quantities; 


<I>  =  sum  of  the  angles  that  adjacent  segments  made  with 

each  other  then  considered  as  vectors  located  in  first 
and  second  quadrant  with  respect  to  a  preferred  direction  y. 

^  =  the  absolute  value  of  the  signed  sum  of  the  angles  that  the 

segments  make  with  a  preferred  direction  when  located  as  above. 


Further  consider  the  conditions  of  the  road  net  in  a  cell.  Because  of 
such  factors  as  a  bridge  out,  road  demolished,  or  mud,  traffic  may  be 
stalled.  From  these  three  quantities  can  be  determined  estimates  of  a 
transfer/no-transfer  function  for  each  cell.  Thus,  as  in  the  case  of 
off-road  traffic,  a  density  p^v2  and  a  mass  absorption  coefficient  k  cv2 
can  be  evaluated. 


Scattering  by  Enemy  Action 

k^yj.and  k£jv2  are  just  the  probabilities  per  unit  occurrence  that  a 
scattering  event  will  take  place. 


Scattering  by  Terrain 

k  ev[  is  the  probability  per  unit  no-go  that  a  scattering  event  will  take 
place.  The  on-road  density  pevz  is,  however,  in  general  different  from 
p bv2  ■  Consider  Figure  9;  a  scattering  event  occurs  if  there  are  two 
or  more  branches  issuing  from  a  cell  with  different  tangent  vectors. 

We  will  assume  that  the  case  in  which  two  road  branches  issue  from  a 
cell  is  fundamentally  the  same  as  that  in  which  multibranches  issue  from 
a  cell.  A  further  refinement  would  be  to  make  this  distinction.  Thus 
we  could,  in  principle,  determine  a  scatter/no-scatter  function  for  each 
cell  and  then  derive  a  density  pev2  and  corresponding  coefficient  k  ev2  . 


GENERAL  SIGNIFICANCE  OF  THE  PHASE  FUNCTIONS 


The  phase  function  P  was  defined  to  be  the  probability  per  unit 

2tt 

— >  — >  — > 

solid  angle  about  12 '  at  R  that  a  unit  absorbed  at  R  will  be  scattered  into 
— > 

the  12  direction  (see  equations  1.2,9  arid  1.2.  10  ) . 


OPERATIONAL  DEFINITIONS  OF  THE  PHASE  FUNCTIONS 

We  now  consider  separately  the  different  phase  functions  given  in  the 
table . 
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((3)  Speed  Transfer  by  Enemy  Action 


Here  the  on-road  transfer  phase  function  differs  from  the  off-road 
transfer  phase  function.  Realistic  plan  functions  can  be  constructed  in 
the  following  manner: 

First  recall  that 


p(  £2 ,  £2'  ) 


2  7T 


doo  -  the  probability  per  unit  solid  angle  about  £2' that 


the  unit  which  was  absorbed  will  be  emitted  into  dui  about  £2  ,  given  that 
an  absorption-scattering  (transfer)  event  has  occurred. 


Now,  consider  pa 
in  F igure  10  . 


(v2  _>Vj  ,  £2,  if')  and  consider  the  situation  illustrated 


Once  a  transferring  event  has  occurred  on  the  road,  a  unit  traveling  in 
the  £2 1  direction  becomes  reoriented  in  the  £2  direction,  with  a  symmetric 
distribution  about  the  preferred  direction  of  travel,  x.  As  with  the  phase 
function  for  pure  scattering  to  be  considered  later,  we  assume  that 
pa  (v2_>  v(  ,  if,  if  1 )  =  pa  (v2_> V[  ,  if),  or  the  phase  function  is  independent 
of  if1  .  Thus, 


Ztt  Pa 

(v2— >Vj 

£2  )  I  v2 

— )  — ) 
(R,  £2 

')  dco' 

Pa 

(VZ-lVj, 

0)  /„' 

hz 

(R,  O')  da>' 

Ztt 
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Now  consider  pa  (vj  ->v,  ,£},  Q')  and  consider  the  situation  illustrated 
in  F igure  1 1 . 

x  Preferred  Direction  of  Travel 


Figure  11. 

— > 

Once  a  transferring  event  has  occurred,  a  unit  traveling  in  the  S21 
direction  becomes  reoriented  in  the  Q  direction  independent  of  fi1.  Let 
the  be  the  tangent  vectors  to  the  roads  as  explained  in  paragraph  e 
under  "S^peed  Transfer  by  Terrain",  where  the  road  has  been  "linearized" 
and  let  fl  be  the  resultant  on-going  direction  of  the  units.  We  assume 
that  the  phase  function  is  of  the  following  form: 

Pa  (vi  i  =  £  5  a  (T-  (R),  SI) 

i  =  i 

where  _ 

<5a  (Tj  (R),  SI)  = 


H  is  the  number  of  road  branches,  and 


-■»  -> 


A  realistic  determination  of  the  oa  (T'  ,  fi)  is  indicated  in  Figure  12. 

x 


Figure  12. 
20 


when  <5  a  (Tj  ,  fl)  is  to  6  a  (T2  ,  il)  as  hf  is  to  h2 
and  where  /u  j  6a  (Tj  ,  ft)  +  6a  (T2  ,  ft)  j  dw  =  1 . 
Thus, 

~  /„•  Pa  (v,  -»  v,  ,  ft,  ft1 )  Ivl  (R,  S2 1  )  dw* 
2.5.2  /u,  ly,  (R,ft')  dw'  . 


2rr 


Here  the  on-road  case  differs  from  the  off-road  case. 

Off  Road:  We  have  determined,  considering  that  a  unit  is  goal- 
oriented,  that  a  realistic  phase  function  can  be  determined  in  the  following 
manner:  To  fix  discussion,  let  the  positive  x  direction  in  a  rectangular 
system  be  the  "ideal"  direction  of  traffic  (see  Figure  13). 


Normal  Curve 

Figure  13. 
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Assume  that  once  a  scattering  "event"  has  taken  place,  the  unit  independent 
of  ft'  seeks  to  orient  itself  in  the  positive  x  direction  resulting  in  an  ft 
orientation. 

Thus, 


->  ->  -> 


therefore  we  obtain 


On  Road  :  To  a  first  approximation,  on-road  scattering  by  enemy 
action  does  not  take  place. 


Therefore, 

2.5.7  fu j'  Pcv2  (ft ,  ft  1 )  I  v2  ( R i  ft 1  )  dwi  =  0 . 

(e)  Scattering  by  Terrain 

Here  the  on-rcad  case  differs  from  the  off-road  case. 


Off  Road:  We  assume  that,  fundamentally,  the  phase  function  for 
terrain  scattering  is  of  the  same  form  as  that  for  enemy  action. 

Thus , 

2.5.8  ^  Ju.  Pdvi  (ft)  Jv.  (R,ft')  dJ=  Pdwjil  /w,  Iyl  (R.ft'Jdco'. 

On  Road :  The  on-road  scattering  phase  function  might  be  viewed 
as  in  the  cases  of  speed  transfer  due  to  enemy  action  and  due  to  terrain. 
Consequently,  we  obtain 

2.5.9  27  Ao'Pdvz  (ft.ft')Iv2  (R,  ft !  )  du>'  = 

M  —  ->  -> 

2  6  (T;  ,Si) 

-  /„.!«  (R.B')du,'. 
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THE  MODIFIED  GENERAL  EQUATIONS 


Clearly,  the  definitions  and  assumptions  introduced  in  this  chapter 
mathematically  simplify  the  general  transfer  equations.  The  appropriate 
simplified  equations  will  be  derived  below. 

Time-Independent  Equations 

The  absorption  terms  given  in  the  section  "Notes  on  the  Specific  Transport 
Equation"  in  Chapter  1  are  unaltered.  The  emission  terms  will  take  on 
the  following  form: 

■Wi  <v!  -vi  )  =  Pa  <v2  ^v]  '^)  d“' 

W  K-v,)  =  ^  S  «a  (T;  ,3)/u,Ivl  (R,fbdu' 


Jpvi  (V2  ->V1  ) 


Pb  (v2  ->vi.  /w'  !V2 


Jp  v2  ( vi  — »  v2  )  =  2t r  ^  b  (^i  >  ^ 1  /co'  ^  vi  )  di°' 


KH  V  1  - >  f  - 7  — 7 

V  =  -?r  pCvi«a>  ^V, 


Jw2  -  ^ 


Pdvi  (fi)  O  1  vi  n')du' 


k  pxrp  M  ^ 


2  6d  (Tx  ,fi)  /  ,  IV2  (R,n')da,'. 


— ?  — > 

d  Iwi  (R,  & )  — >— >  .  — >  — > 

— ] ^nr - '=  -  Kvi  ivi  (R.fl)  +  stvi  /«'  ivz  (R.«')d“' 


+  ^SVl  / U) '  I  VI  (R>  ^  ')  du>' 


AND 


d  I  V2  (Ri  fi) 


=  -K^  Iw  (R.fl)  +  Stv2  /ui  IV1  (R.Q)du)' 


+  2sv2  /u.  Iv2  (R.O)d«' 


where 


-  ^bv2  Pav2  /  7?  >  ,  ^  cv2  Pbv2  , 

"  - 2F~ - Pa  (v2  -*  v,  Q  )  +  - ^ -  pb  (v2  _>  vt 


kdvi  Pavi  kevi  Pbvi  .7*. 

Pcvi  (A)  +  2-n-  Pdvi 

^-P-a-V1-  Z  «a(Ti,0)  +  Z  6b  (Tj  ,  0) 

2tt  i  =  l  2tt  i  =  i  D 

kevz  Pevz  £  r  F. 

— S -  ,5,6d(Ti  .a)- 


Time-Dependent  Equations 


The  time -dependent  equations  become: 


dlvl  (R.fl)  _  dlvl  (R,  S2 )  v  T  7l\ 
dt  d^  "  Kvl  Ivl  (R’n) 


d  I  v2  ( R  i  fl ) 


+  Stv.  L'  1 V2  W.B'l-W  +Ssv,  /„'IV1  (Hi  a1  )du' 

■  -  _kv2  ivz  (3.3) 

as 


+  2tv2  /u'  I  vl  (R,  fi )  dco  +Sgv2  /w  i  Iv2  (R,fl  )  du> 
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la  general,  we  will  be  interested  in  the  density  functions  that  are  conse¬ 
quently  governed  by  the  following  equations; 


2.  6.  3 


where 

u,  =  u  ,  (R)  u2  =  u2  (R) 

2n  =  —  kyz  Pavz  +  ^  cv2  Pbvz 

2ir  2tt 

v  _  ^dvi  Pavi  i  ^evi  Pbvi 

^S1  "  — zli -  - Zi - 


where  arbitrary  source  and  loss  terms  have  been  dropped. 
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